Late time solutions for inhomogeneous ACDM cosmology, their characterization and 

observation 
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Assuming homogeneous isotropic ACDM cosmology allows A, spatial curvature and dark matter 
density to be inferred from large scale structure observations such as supernovae. The purpose of this 
paper is to extend this to allow observations to measure or constrain inhomogeneity and anisotropy. 
We obtain the general inhomogeneous anisotropic ACDM solution which is locally asymptotic to 
an expanding de Sitter solution as a late time expansion using Starobinsky's method (analogous 
to the 'holographic renormalization' technique in AdS/CFT) together with a resummation of the 
expansion. The dark matter is modeled as perfect dust fluid. The terms in the expansion system- 
atically describe inhomogeneous and anisotropic deformations of an expanding FLRW solution, and 
are given as a spatial derivative expansion in terms of data characterizing the solution - a 3-metric 
and a perturbation of that 3-metric. Leading terms describe inhomogeneity and anisotropy on the 
scale set by the cosmological constant, approximately the horizon scale today. Higher terms in the 
expansion describe shorter scale variations. We compute the luminosity distance-redshift relation 
and argue that comparison with current and future observation would allow a partial reconstruction 
of the characterizing data. We also comment on smoothing these solutions noting that geometric 
flows (such as Ricci flow) applied to the characterizing data provide a canonical averaging method. 



I. INTRODUCTION 

Modern large scale structure observations, such as the 
SDSS low redshift studies of 'standard candle' supernovae 
[l|, |2|, afford us the capability to directly test our cos- 
mological model assuming only the Einstein equations, 
matter content relevant for our current epoch, and ho- 
mogeneity and isotropy (and hence the FLRW metric). 
Historically Hubble's measurement of the current expan- 
sion rate was the first test of this type. Recently the 
famous results of the Supernova Search Team 3] and the 
Supernova Cosmology Project [4] allowed comparison of 
luminosity distance-redshift for various matter content 
models with supernovae data up to redshift ~ 5, showing 
the likely current domination of a positive cosmological 
constant over dark matter. 

An important point is that unlike CMB observations, 
large scale structure allows inferences to be made about 
the parameters and dynamics of a model without any 
assumption on initial conditions, except for the, albeit 
strong, assumptions of homogeneity and isotropy. While 
the cosmological constant can be estimated from CMB 
measurements, this inference is far less direct, and is sub- 
ject to a choice of inflationary initial conditions. Here we 
are certainly not arguing that one should not take infla- 
tionary initial conditions seriously, but rather we wish to 
emphasize that when one makes inferences about physi- 
cal models from observational data, the strength of the 
parameter estimation implicitly depends on the strength 
of the assumptions made. It is useful to have a variety of 
methods to allow estimation with differing assumptions. 

It is worth emphasizing that in the context of infla- 
tion it is the inflationary mechanism that suppresses large 
scale inhomogeneity, anisotropy as well as spatial curva- 
ture. Thus if one is interested in assessing the magni- 
tude of spatial curvature by measurement from observa- 



tion, for example if one thought inflation lasted only a 
minimal time and therefore interesting structure may be 
seen on large scales (e.g. [^]), it is conceptually incon- 
sistent to include spatial curvature. One should consider 
general anisotropic inhomogeneous deformations at the 
same time as introducing spatial curvature, unless one 
supplements the initial conditions before inflation with 
homogeneity and isotropy, as for example occurs after 
semiclassical bubble nucleation of an open universe [6| . 

The purpose of this paper is to allow the comparison of 
a ACDM cosmology with large scale structure data such 
as supernovae, without making the assumptions of homo- 
geneity or isotropy on large scales. This is not because 
we are arguing the universe is far from homogeneous or 
isotropic, but simply because in order to argue the uni- 
verse is homogeneous and isotropic, ideally one should 
deduce this from the data without assuming it. Fur- 
thermore with future measurements one would hope to 
be able to quantify precisely how far from homogeneity 
and isotropy the universe is, or at least put rigorous up- 
per bounds on how large these quantities are on different 
scales. 

Of course in order to do any of this, one must be able 
to compute and characterize general inhomogeneous solu- 
tions for a ACDM cosmology. In addition, one must then 
be able to compute quantities relevant for observation, 
such as the luminosity distance-redshift relation, and 
show how this allows inferences about the data that char- 
acterizes the solutions to be made from real large scale 
structure data. Certain exact solutions have allowed such 
analysis to be performed for very special cases, for in- 
stance using the LTB solution which is isotropic but inho- 
mogeneous, or the Bianchi models which are anisotropic 
but homogeneous |7Hl3|. but general exact solutions are 
likely not to exist in any closed form. 

In this paper we show how to construct general inho- 



mogeneous, anisotropic solutions of the Einstein equa- 
tions with cosmological constant and dust fluid as a 
model for baryonic and dark matter using a method 
which we call a 'late time expansion'. This method was 
first developed by Starobinsky jlj| in the context of in- 
fiation, and is similar to methods eniployed in AdS/CFT 
holographic renormalization (see fl»|). Rather than 
starting from the conventional approach of cosmological 
perturbation theory about the FLRW solution, instead 
we use the fact that provided the local universe is dy- 
namically dominated by the cosmological constant, one 
can describe the solution in a precise way as a deforma- 
tion about de Sitter space, characterized by its late time 
asymptotic behaviour, or what we shall term the final 
boundary conditions. The leading late time behaviour of 
the solution is a universe that is locally de Sitter. Choos- 
ing the late time de Sitter conformal boundary to have 
proper time, i — >■ oo, the full solution then takes the 
form of an expansion in positive powers of e~ '^, where 
higher terms in the expansion give inhomogeneous cor- 
rections to this local de Sitter behaviour, and are given 
in terms of increasing numbers of derivatives acting on 
the final boundary data. We find the form of this final 
data is given by a 3-metric and a perturbation of that 
3-metric. The leading terms in the late time expansion 
describe inhomogeneity and anisotropy on the scale set 
by the cosmological constant, and hence for our universe 
would correspond to approximately today's horizon size 
deformations. Higher terms in the series then govern 
derivative corrections which correspond to increasingly 
smaller scale inhomogeneity and anisotropy. 

The first part of our paper is devoted to giving the so- 
lution for ACDM, taking dust fiuid, using the 'late time 
expansion' about an asymptotically locally de Sitter so- 
lution with fiat spatial slices. We give results beyond 
leading order, and in particular discuss the final data 
that arises in the expansion and characterizes the solu- 
tion. For dust fiuid it is possible to take a normal coor- 
dinate system adapted to the conformal boundary that 
is comoving with the dust - this is not possible for fiuids 
with pressure. This simplifies the calculation consider- 
ably, and is what allows us to give results straightfor- 
wardly beyond the leading order. 

In the second part of the paper, having showed how to 
give inhomogeneous anisotropic solutions as a late time 
expansion deforming about fiat sliced de Sitter, we go 
on to consider the convergence of the late time expan- 
sion, arguing that whilst one might naively expect it to 
break down around A-matter equality, it actually is far 
better, extending back to very early epochs. We then 
discuss the drawback of our first expansion, namely that 
homogeneous isotropic dust and curvature are also de- 
scribed as deformations and treated order by order in 
the late time expansion, even though we know the exact 
homogeneous isotropic solution is an expanding FRLW 
universe, whose solution we may treat analytically. How- 
ever, we show that the result may simply be modified to 
give an expansion about the exact expanding FRLW so- 



lution using a simple but powerful series resummation 
on the late time expansion. The terms in the expansion 
then describe deformations only involving inhomogeneity 
and isotropy. If the deviations about FLRW solution are 
small we believe, although here do not prove, that this 
late time expansion converges back to high redshifts, cer- 
tainly sufficient to describe the region relevant for large 
scale structure observations. 

The third part of the paper is then concerned with ap- 
plying these inhomogeneous and anisotropic solutions to 
observations of our universe. We compute the luminosity 
distance-redshift relation explicitly for the first two non- 
trivial orders in the expansion. We then discuss how one 
might in principle use this to compare to supernovae data 
and therefore extract information about the characteriz- 
ing data, the boundary 3-metric and its perturbation. In 
particular we discuss how the boundary data is extracted 
as a derivative expansion about the point corresponding 
to the observers asymptotic position, and that fitting in- 
creasingly high orders in the luminosity distance-redshift 
relation expansion to data gives information about in- 
creasingly high derivatives of the local curvature of the 3- 
metric and perturbation. Thus information characteriz- 
ing the inhomogeneous and anisotropic universe we live in 
can be systematically extracted from observational data, 
provided that, as we argue, we are living in and observ- 
ing the regime where the late time expansion converges. 
In fact not all characterizing data can be extracted from 
this luminosity distance-redshift measurement, and we 
deduce the subset of the data that may be extracted. 

We briefly conclude with a discussion of averaging or 
smoothing of our inhomogeneous solutions. Since in our 
method the late time inhomogeneous cosmologies are de- 
termined by geometric final data consisting of a 3-metric 
and its perturbation, geometric flows, such as Ricci flow, 
provide the canonical notion of averaging this character- 
izing data. Considering this we arrive at a natural av- 
eraging of our solutions similar in spirit although rather 
different in origin to that discussed by Buchert and Car- 
fora [la, [TtI • It is simple to argue that provided one be- 
gins with small deformations of FLRW, such smoothing 
flows to the FLRW solution. 



II. LATE TIME EXPANSION ABOUT DE 
SITTER 



Following the method of Starobinsky's late time ex- 
pansion |14| , utilized in the context of inflation, we now 
show how to construct the general cosmological solution 
in a ACDM model that is locally asymptotic to de Sit- 
ter, taking the dark matter and baryonic matter to be 
modeled by a single dust fluid. 

Unlike a conventional treatment of inhomogeneity as a 
deformation about de Sitter in a perturbative expansion 
in deformation amplitude, here Starobinsky's method 
treats the solution as a deformation about de Sitter as a 
late time expansion, where the perturbation parameter 



instead measures how late the region being described is. 
No assumption of homogeneity is necessary, and inhomo- 
geneity is treated fuhy non-perturbatively. 

A key technical tool that we may use is to choose a con- 
formal frame where the dust is comoving. In fact this is 
possible only for dust fluid, and very much simplifies the 
construction of the solution, which previously was given 
by Starobinsky and later more rigourously by Rendall 
[18j for cosmological constant and general perfect fluid 
but only to leading order. Here, for dust, we are able to 
straightforwardly derive several orders in the expansion. 

We proceed by writing the metric covering the late 
time region of interest in a Gaussian normal form to con- 
stant y surfaces. 
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where, as before, a;' for i = 1,2,3 give spatial coordinates 
on the constant y surface and H^ = A/3, and would give 
the Hubble rate for the de Sitter geometry corresponding 
to the cosmological constant A. By construction of nor- 
mal coordinates the vector d/dy is tangent to timelike 
geodesies normal to the constant y surfaces. We have 
the freedom to choose these constant y surfaces to co- 
move with the dust fluid so that the stress tensor takes 
the simple form, 

'^yy = Jp:^ (^ + p) - Tyi = o, T^ = -^.9y (2) 

and this fixes our choice of conformal frame. Taking a 
different frame would require us to introduce a velocity 
field for the dust fluid. 

We decompose the Einstein equations (taking SttG = 
1) to give a scalar and vector equation with respect to 
the 3-metric gij, 

V^g^J-V^g = 
and in addition a tensor equation, 
2 1 
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where indices are raised and lowered with respect to gij, 
Vi is the Levi-Civita connection of gij, Rij is its Ricci 
tensor, and we define gij = dygij, g — g^-' gij and likewise 
cjij = dygij, g = g^-'gij- As usual, the fluid dust equation 
of motion is implied by the Einstein equations due to the 
contracted Bianchi identity. 

We now consider inhomogeneous solutions and use 
a late time expansion as a tool to solve the equation. 
Firstly we consider solving the equation ^ as an expan- 
sion in y. One obtains, 

9^J {y,x) = gij {x) -f y'^af^ {x) + y^hij (x) 

+y*a^^{x)+y'a^\x) + ... (5) 



where the two pieces of data for the solution of this sec- 
ond order equation are provided by the symmetric gij (x) 
and hij{x). We may interpret gij{x) as the metric of a 
3-geometry - the conformal boundary metric of the full 
spacetime - and then the other data hij (x) we regard as 
a symmetric tensor living on this geometry. The terms 
al,j and all higher terms in the expansion are determined 
in terms of this metric and the data hij . Explicitly we 
find these lower order terms to be, 
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where Rij is the Ricci tensor of the data gij and indices 

(i) 

in the terms a^ are lowered and raised with respect to 
gij, and Vi is its Levi-Civita connection, and O is a dif- 
ferential operator that we give below in equation (|IOI) . 

We now turn to the remaining scalar and vector equa- 
tions. The scalar equation simply determines the be- 
haviour of the dust fluid density as. 
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^3y^h + ya^''^h + 0{y'^) (7) 



The vector equation is non-trivial and constrains the 
space of solutions. Let us define the vector $i = V^g^ — 
Vij, and then assuming the tensor equation Q is satis- 
fied we find, dy^i = ( ~ ^ 55) *^i which may be solved 
for each component of <i>i separately by quadrature to 
give. 



$,(y,a;) = Aj(x) y2e-i/o'''J'ste.^) 



(8) 



with integration constants Ai{x). However, evaluating 
$i on our late time expansion above, we find. 



My, a;) = 3 (V^7iy - W^h) y^ + 0{y^) 



(9) 



and hence we obtain Ai{x) = 3 (S/^hij ~ \7ihy Thus 
in order to set the vector equation to zero, we require 
that Ai(x) = 0, and hence we require a constraint on our 
initial data hij, V^hij — Vih = 0. 

The form of the operator O, which we previously post- 
poned giving, is. 



o,^"'"T^n - --v2r,, + -v,9,r-i?™"r„ 



+Rm(iTj)'^ - -tRT^j + -gij RrnnT'"'^ (10) 

for action on a symmetric tensor Ty which obeys 'S/^Tij — 
ViT = as do both «„„ and /im„. We note the fact 



that in 3-dimensions the Riemann tensor which enters 
the above equation is simply determined by the Ricci 
tensor as, 
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due to the vanishing of the Weyl tensor. 

This 'late time expansion' method was first used by 
Starobinsky in the context of inflation to argue that 
inhomogeneity is washed out during inflation. In |14l | 
Starobinsky stated the result for cosmological constant 
and general fluid matter, but only to leading order. Re- 
cently Rendall [18] has given a more rigorous derivation 
for Starobinsky's somewhat briefly justified leading or- 
der results. To our knowledge, Starobinsky's work has 
not been extended beyond leading order for fluids, and 
has never been applied to the late time universe, but only 
to inflation. 

We note that Starobinsky's method is quite analogous 
to the techniques developed in geometry for hyperbolic 
space [19] which are ubiquitous in string theory hologra- 
phy (for a recent review see [3). Note that in [l^ [2^ 
beyond leading order results are given for pure cosmo- 
logical constant matter. There has been considerable in- 
terest in using string theory inspired holographic ideas 
to give new quantum gravitational theories of the early 
universe and inflation [21-24]. We emphasize that whilst 
some technical methods of our current work bear similar- 
ity with those used in holography conceptually our work 
is very different in character. We are here concerned only 
with describing the late time universe and not its early 
stages, and have no dual holographic description in mind. 

We also note that the late time expansion method may 
apply more generally than to cosmologies that are asymp- 
totically locally de Sitter. For example, in [23 analogous 
expansions have been proposed for asymptotics with fu- 
ture singularities. One might hope that one can also 
generalize this method to allow a late time expansion for 
a cosmology with a quintessence field and dust, and this 
would be an interesting direction for future work. 



A. Characterizing solutions 



for A € K"*". Defining the equivalence, hij ~ A^/iy for 
A G M+, we may then take the data characterizing the 
solution to be given by the boundary metric and equiva- 
lence class of hij, namely (gij, [hij]). 

Consider a 3-geometry with metric gij and an infinites- 
imally perturbed geometry, gij+eSgij, for infinitesimal e. 
Due to the ability to rescale the infinitesimal e, we take 
a geometry and infinitesimal perturbation to be charac- 
terized by the pair {(jij, [Sgij]). Define the vector field, 
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(14) 



Consider a coordinate system such that the 'harmonic 
gauge' condition Ui = holds. This locally fixes all the 
coordinate freedom. We do not address whether there 
are global obstructions to the existence of such a gauge, 
but believe it is plausible that such a gauge always exists 
and is unique, and assume so from now on.[28J In such 
a gauge, taking hij = Scjij — \gij6g, then hij obeys pre- 
cisely the appropriate condition ([12]). The physical dust 
density constraint, h < 0, then becomes 6g < 0. Since 
S ( -^det (jij ) = 2 -v/det^ijJ^, this implies the perturba- 
tion in harmonic gauge locally decreases 3-volume. 

Thus we conclude that a 3-geometry and perturbation, 
{(jij, [Sgij]), can be mapped to the data (gij, [hij]) of our 
late time solution, provided the perturbation can be cast 
in harmonic gauge Ui = 0. A physical perturbation in 
that gauge locally decreases volume, which ensures posi- 
tive dust density. This 3-geometry and perturbation pro- 
vide an elegant geometric way to characterize our late 
time inhomogeneous solutions. 

We note that if one wished to include other matter in 
addition to the dust fluid, then one would have to intro- 
duce additional tensor fields to describe its data. As we 
know from discussion of holographic renormalization in 
AdS/CFT, extension to include other matter fields would 
likely be a straightforward, but messy exercise. We will 
not be concerned with that here, as from a practical point 
of view we wish to use this solution to better understand 
our late time universe, and other forms of matter (radi- 
ation, infiaton etc) play little role at late times. 



Our choice of conformal frame makes it straightforward 
to identify that our late time solution is given by the 
expansion ([S|) where the data is the (Euclidean signature) 
3-geometry, specified by the metric gij, together with a 
symmetric tensor hij defined on this 3-geometry obeying. 



y'h. 



Vih^O . 



(12) 



The late time density of the dust is given in terms of 
the trace, p{y,x) = —3Hh{x)y^ + 0{y^), and hence we 
require this trace to be negative definite h = g^-'hij < 
to ensure a physical dust density. We note that there is 
a global scaling invariancc of the metric ([TJ , 



y~^y/X, x'-^xyX, g^^ 



h,j ^ X^ h,j{lZ) 



B. Remarks on the late time expansion 

We have seen that the terms in our late time ex- 
pansion give a derivative expansion in the boundary 
data {gij, hij). The derivative expansion is controlled by 
the cosmological constant scale, 1/H . Hence the lead- 
ing terms in the expansion describe inhomogeneity and 
anisotropy on the largest scales, namely the late time 
horizon size. Higher terms in the expansion then describe 
shorter scales, and in principle the derivative expansion 
may describe all scales, although this is contingent on the 
late time expansion converging. In general, one might 
expect that as one moves to earlier times so that the in- 
homogeneity becomes significant on scales smaller than 



the comoving horizon size one would expect convergence 
to worsen and at some point break down. 

The expansion need not describe the global future of 
the spacetime, but may just describe a local region of the 
future conformal boundary that is asymptotic to de Sit- 
ter. Other behaviours are possible in the future, such as 
dust or curvature collapse, corresponding to a future con- 
formal boundary that is a black hole singularity. Then 
the hyperbolic nature of the Einstein equations presum- 
ably implies that the late time expansion converges only 
within some portion of the past light cone of a region that 
is asymptotic to de Sitter in the future, as illustrated in 
figure [T] However, we note that for our universe, on large 
scales, we do expect that the universe is asymptotic to de 
Sitter, provided that the matter content is a cosmologi- 
cal constant, and dark and baryonic matter is effectively 
describable by dust fluid. While locally matter may col- 
lapse to form black holes, one should average over such 
small scale inhomogeneity when discussing inhomogene- 
ity with size comparable to today's horizon size as our 
late time expansion does. 

One might wonder what the essential difference is be- 
tween our solution and a Taylor expansion of data about 
any Cauchy surface. In effect one can view our solution 
as a Frobenius expansion about a regular singular point 
in the Einstein equations, that given by y — >■ 0. It is 
worth emphasizing the following points. 

1. Whilst one can indeed characterize a solution by 
analytic data on any Cauchy surface and then give 
the analytic expansion governing its evolution (in 
the sense of Cauchy-Kowaleski) there is no canon- 
ical surface to perform this about. In our con- 
struction the future conformal boundary provides 
a canonical surface. In particular this means that 
through making observations one can in principle 
measure model parameters which are independent 
of the observer. 

2. Data to the future of a region of a Cauchy surface 
is not guaranteed to enter a cosmological constant 
dominated future - for example, dust or curvature 
may come to dominate and lead to a singularity, as 
illustrated in figure [1] By characterizing our solu- 
tion with data that lives at the future de Sitter con- 
formal boundary, we have focused by construction 
precisely only on the regions of cosmology where 
the cosmological constant comes to dominant, as 
we as observers see today. 

3. On a more technical level, the late time expansion 
exhibits a separation of the orders at which var- 
ious matter/curvature contributions first appear. 
For example, the spatial curvature of the bound- 
ary metric first enters at order y^ . Dust first enters 
the expansion at order y'^ . Had we included ra- 
diation, it would not enter until order j/^. This is 
generally not the case for a Taylor expansion about 
some Cauchy surface, where all forms of matter and 



curvature enter at the first order in the expansion. 
We believe that this clean initial separation at low 
orders in y will prevent possible ambiguity when 
fitting data, particularly at low redshifts. 

We now bricfiy discuss our expansion from the per- 
spective of AdS/CFT holographic renormalization. If 
a dual QFT existed, it would live on a 3-d spacetime 
with metric gij. The tensor hij would be interpreted 
as giving the boundary theory conserved stress tensor 
as Tij = hij — (jijh. Due to the presence of the dust 
matter, the stress tensor is not trace free, and hence 
the dual theory is not conformal. The dust would be 
associated to a marginal dimension 3 operator in the 
dual theory, and a trace perturbation of gij and hij give 
the non-normalizable and normalizable gravity modes re- 
spectively that correspond to the source and vev of this 
operator. An interesting point is that there are no ir- 
relevant operators, only the marginal graviton and dust. 
The absence of irrelevant operators corresponds in the 
gravity to our late time expansion being a stable attrac- 
tor solution with no perturbative unstable modes that 
might destroy the de Sitter late time asymptotics. 



C. Convergence and application to our universe 

Using the late time expansion method we have shown 
how to construct the general inhomogeneous anisotropic 
late time solution as a deformation of de Sitter. Then 
one expects that provided one lives at late times, so that 
the cosmological constant dominates, then the expansion 
would converge, and one could truncate the expansion to 
obtain a controlled approximation to the true solution. 
From the practical perspective of application to our uni- 
verse a natural question follows. Do we live late enough 
for the series to converge? Recall that in our expansion 
above the homogeneous dust component is also treated 
as part of the deformation about de Sitter. Whilst we 
expect A to dominate matter today (as flm ~ 0.3), it is 
not parametrically larger. Naively one might worry even 
whether the series converges restricting to FLRW for us 
today, or certainly back to earlier redshifts, say z ^ 5, 
where supernovae yield important large scale structure 
observations. At such an epoch matter dominates A for 
FLRW and thus one would naively assume this is too 
early an epoch for our late time expansion to describe. 

This reasoning however is too quick as we shall now 
show. Rather surprisingly, for FLRW relevant for our 
universe, the expansion converges back to the big bang 
(in this dust- A model). Let us consider a fiat sliced 
FLRW universe. The solution is simply written using 
our form of the metric as, 
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FIG. 1: A cartoon of the regions where the late time expansion might converge for a cosmology with both de Sitter and black 
hole future conformal boundary asymptotics. 



form of the exact solution that convergence of the se- 
ries in fact will extend back to the big bang time, as 
this represents the only non-analytic behaviour (a branch 
point) in the solution. To approximate the solution well 
at early times, of course one will require many terms to 
reach a given precision, compared to approximation at 
late times. However, provided one takes enough terms, 
the series will sum to give the required level of approx- 
imation. Thus we see, at least for flat sliced FLRW, 
that the late time expansion actually works rather better 
than expected in the sense that it converges for all the 
history of the universe, back to the epoch where other 
matter forms become relevant. Adding spatial curvature 
to the FLRW model appears not to change this picture, 
although we have not proved this, but have observed in 
examples that convergence extends back to the big bang, 
even for large curvature components. 

We now make a brief aside on radiation which we have 
not included in our model. Radiation enters the late time 
expansion at order y^. In principle we should include 
this when discussing terms in the expansion of order y* 
or higher. However, we note that since matter-radiation 
equality is at such an early epoch, z ~ 3000, the numeri- 
cal size of the terms controlling radiation in the expansion 
must be unnaturally small [29i], of order ^ 10^^^, com- 
pared to the size of the terms controlling matter (such 
as po in the exact FLRW solution above). Thus whilst 
radiation in principle enters the expansion at order y*, 
the size of the contributions would be very suppressed 
compared to the size of terms for matter and curvature 
provided one is studying late times in the universe, i.e. 
much after matter-radiation equality. 

We have arrived at the rather remarkable conclusion 
that provided inhomogeneity and anisotropy are small, 
then our late time expansion presumably converges all 
the way back to the beginning of the matter era. From 
the perspective of using the late time expansion to model 
observations this is very encouraging, as presumably tak- 
ing sufficiently many terms, one can obtain accurate ap- 



proximations to the true solution far enough back to 
include all current (and future) supernovae data, even 
though high redshift data fall outside the cosmological 
constant dominated epoch. 

However, having found that our solution converges 
much further back than A-matter equality we then see 
that actually the fact that we have deformed about the 
late time de Sitter solution is somewhat unattractive, in 
the sense that we only describe homogeneous dust ap- 
proximately. One might expect that if inhomogeneity 
and anisotropy were small one would require only a few 
leading terms to accurately model our universe. However 
looking back to the matter epoch our late time expansion 
will require many terms in order to accurately approxi- 
mate FRLW. 

We may phrase this another way. Corrections in inho- 
mogeneity and anisotropy enter at the same level as the 
homogeneous dust component. Whilst we have separated 
the dynamical scale of A from those of inhomogeneity and 
anisotropy, we have not managed to also separate out ho- 
mogeneous dust. We may see this concretely exhibited in 
the solutions above for the expansion ([S|) as anisotropic 
terms (those tensors not proportional to Sij), and inho- 
mogeneous terms (those involving derivatives) enter the 
expansion in the same way as isotropic homogeneous spa- 
tial curvature contributions, which go as R6ij, and dust 
terms, going as hSij. 

This situation may be simply improved using the fact 
that since we know the homogeneous isotropic solution is 
simply the FLRW solution, we may directly factor it out 
of the expansion, leaving a late time expansion that only 
consists of terms that explicitly involve inhomogeneity or 
anisotropy. We may say that instead of deforming about 
late time de Sitter, we instead rewrite our solution as a 
deformation about late time expanding ACDM FLRW. 
In the next section we proceed to implement this. 



III. IMPROVEMENT BY RESUMMATION: A 
LATE TIME EXPANSION ABOUT FLRW 

Consider a general expanding FLRW solution writ- 
ten using our above late time expansion. Then we have 
gij = a^{y)Q.ij, where fly is the round sphere metric for 
a closed universe, the flat metric 5ij for flat spatial slices, 
or a hyperbolic metric for an open universe. In terms 
of the characterizing data, gij — Qij with Rij — 2k gij , 
with fc > 0, fc = or /c < for the closed, flat and 
open cases respectively. Let us take the constant po to 
characterize the dust density relative to the cosmological 
constant and spatial curvature. From above we see the 
dust is characterized by the trace of hij. Since we are 
interested in a homogeneous and isotropic solution, then 
we see that hij takes the form, hij = — p<^ / (2>H)'^Q,ij , so 
that asymptotically p ^ poy'^ + . . .. Thus wc find that 
the FLRW solution is given by the late time expansion. 
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By construction we obtain the general FLRW solution, 
with its homogeneous and isotropic spatial curvature and 
dust as a late time series expansion about local flat sliced 
de Sitter. 

However whereas for inhoniogeneous anisotropic solu- 
tions we have no analytic control over the solution beyond 
our late time expansion, for the homogeneous isotropic 
solution we know the solution, simply by solving the 
FLRW equations. Let us denote the full FLRW solution 
above as, gij{y) = a'^{y;k,po)nij, so that the function 
a{y;k,po) is the solution of the equations. 
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with the appropriate asymptotic behaviour p ^ p^y^ -f 
.... The expanding solution to these odes may be 
straightforwardly computed numerically, or given as a 
Frobenius series expansion analogous to (llSp (expression 
(jlSp is without curvature) to very high order. Alterna- 
tively one may manipulate the formal solution written 
implicitly in terms of the integral. 
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expansion did not vanish in the homogeneous isotropic 
limit, namely those that describe the correction to de Sit- 
ter represented by the expanding FRLW solutions. How- 
ever, given that we know the form of the FLRW solutions, 
we can factor out — or 'resum' — precisely these terms 
by taking the metric. 



(18) 



where y = yo gives the position of the big bang. In 
practice the numerical approach is probably the easiest 
method to obtain the function a{y; k, po) to a desired 
precision in terms of the variables y,k, po- 

In our original discussion we presented our late time ex- 
pansion as a power series solution for the metric gij{x, y) 
in equation ([T}. As we discussed above, some terms in the 
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and now performing a late time expansion on gij(x,y), 



mj {Vix) = gij {x) + y'^bfj' {x) + y%j (x) 



(20) 



where we denote the anisotropic part of a tensor Tij (with 
respect to the 3-metric gij) as. 
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1- rj. 



(21) 



Again, the characterizing data is given by the same 3- 
metric gij, and tensor hij that lives on the 3-geometry 
defined by gij, and we note that as before, the tensor hij 
must satisfy the same condition. 



V^/l,; 



V,/i = 



(22) 



and hence can be thought of as being defined by a general 
perturbation of the 3-geometry (subject only to the trace 
inequality for physical dust density) . Now the coefficients 
entering the expansion are found to be. 
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and we see that having resummed the homogeneous 
isotropic terms, we are left only with terms that vanish 
in the homogeneous isotropic limit. 

Using our resummed late time expansion we may ex- 
actly describe expanding homogeneous isotropic ACDM 
cosmologies. Inhomogeneous anisotropic cosmologies 
asymptotic to expanding FLRW ones can be described 
as deformations about FLRW. Since the expansion is 
made in gij , when we approximate a solution by truncat- 
ing the expansion to a given order, we still maintain an 
exact treatment of the homogeneous isotropic dust and 
curvature. We have succeeded in separating the homo- 
geneous isotropic dynamics from the inhomogeneous and 
anisotropic behaviour, allowing the latter to be system- 
atically included as a correction. Since in our universe we 
might expect inhomogeneity and anisotropy to be rather 



small corrections, then one would expect fast convergence 
to the true solution taking few terms in g^ . Then a low 
order truncation of the expansion of gij should give very 
good accuracy back to epochs sufficiently early to include 
current and future supernovae data (say 2 ~ 5). 

We now wish to emphasize an important point. In- 
homogeneity and anisotropy can be treated in cosmolog- 
ical perturbation theory as an expansion in the ampli- 
tude of the metric and matter field perturbation about 
the homogenous isotropic FRLW solution. Our late 
time expansion method also provides an inhomogeneous 
anisotropic solution as a deformation about an expanding 
FRLW solution. The difference of course is that our ex- 
pansion is non-perturbative in the amplitude of the met- 
ric and matter perturbations, but is instead arranged as 
a derivative expansion. Thus, as stated above, we cannot 
treat short wavelengths unlike cosmological perturbation 
theory, without going to high order (providing the scries 
converges). However, we do treat deformations of FLRW 
non-perturbatively in the amplitude of the metric and 
matter field perturbations about FLRW. 

A simple example illustrating this is to take a bound- 
ary metric that is a squashed 3-sphere. Taking a round 
3-sphere would yield spatial curvature, whose effect may 
be small compared to the dust and A today if the cur- 
vature y'^^^R/H^ ^ 1 where we are living at y = j/now 
Obviously, choosing hij appropriately this simply yields 
a closed FLRW solution. However, we may deform 
away from FRLW by taking instead a squashed 3-sphere. 
Again, the effect of the spatial curvature may be made 
small today compared to dust and A. For weak squash- 
ing, the solution can be described by cosmological pertur- 
bation theory about a closed FLRW solution. However 
for strong squashing away from the round metric, the so- 
lution can no longer be described by perturbation theory, 
whereas it can be described by our expansion, and since 
the spatial curvature is small relative to dust and A, we 
would expect our series solution to be convergent back 
to early times in the matter era. 



cone of a late time observer and the redshift of sources 
comoving with the dust fluid. We use this to compute 
the luminosity distance as a function of observer position 
and observation direction, and demonstrate how such 
a measurement encodes the solution data {gij, hij). 
Details of the calculations presented in this section can 
be found in the Appendix. 

A. Late time expansion about de Sitter 



Consider an observer and a monochromatic isotropic 
source of radiation both comoving with the dust. We 
choose coordinates such that the observer is situated at 
(y,x*) = (2/0,0) and the source at {y,x^) = {ye,x],). Now 
consider a null geodesic connecting the two, such that its 
spatial tangent at the observer position is given by u', 
normalized such that v^v^gij{0) = 1. 



The redshift, Z, of a photon emitted at y^ and observed 
at yo may be calculated using the late time expansion in 
powers of y as, 
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(24) 



IV. OBSERVATIONS AND THE 

LUMINOSITY-DISTANCE REDSHIFT 

RELATION 



where V" 



yra;^n _ i^™"(o). The data hj, gij, 



its curvature Rij and higher derivatives are located at 
(y, X*) = (0, 0), i.e. at the observers comoving coordinate 
We now provide a physical example of the application position, but expressed at the future conformal bound- 
of our late time solution. We compute the past light ary. For the luminosity distance, D^, we obtain 
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(25) 



The leading term gives the result for de Sitter with fiat spatial sections. The first non-trivial term then arises 



Future conformal boundary 



h, 




y = ye 



source worldline 

x'{y) = xl 



observer worldline 
xH,y) = 



FIG. 2: The setup for computing the inhomogeneous luminosity distance 



from spatial curvature of the boundary metric. For ho- 
mogeneous de Sitter with closed or open spatial sections 
this yields an isotropic contribution. However, we see 
that for our general inhomogeneous late time cosmology, 
the leading non-trivial contribution depends on the di- 
rection of observation, parametrized by the vector u*. At 
next order there is an isotropic contribution due to the 
trace of hij from the dust fluid, an anisotropic contri- 
bution from hij which we can interpret as being due to 
gravity waves, and further contributions due to inhomo- 
geneity of the boundary metric which involve derivatives 
of the boundary curvature. 

We note that our result here can be interpreted as 
an analog to the luminosity distance-redshift relation 
computed by Kristian and Sachs |26| . where our result 
above is given for a Frobenius expansion about y — 
(our late time expansion), and the result in [26!] can be 
thought of in terms of the local data at the position of 
the observer, which in turn can be thought of as a Taylor 
series expansion about some Cauchy surface containing 
the observers position. Consequently the calculations of 
[26| amount to an expansion in small redshifts. Here, 
our results are valid to all orders in the redshift. As a 
check of our luminosity distance calculation (|25|) in the 
low redshift regime, we have computed the luminosity 
distance using the geometrical optics method of [2g| for 
our late time solutions. We find precise agreement to 
the order of approximation in redshift considered in [26| . 



As a further check we have verified that our expressions 
yield the correct isotropic and homogeneous results in 
the isotropic and homogeneous limit. 



Resummation: Late time expansion about 
FLRW 



The redshift (p4)) and luminosity distance (|25p were 
obtained above using the late time expansion technique, 
where the expansion was performed about de Sitter 
space. We now use the resummation technique described 
in section IHII to improve these results. The key point 
is that even when the homogeneous and isotropic limit 
is taken, both still evaluate to an infinite series. This is 
because we have performed a late time expansion about 
flat de Sitter space. However we know how to obtain the 
solution in this limit, it is just the FLRW solution in the 
presence of dust, curvature and A. We can therefore fac- 
torize our results into a part which gives the redshift and 
luminosity distance for the FLRW solution, and a correc- 
tion written as a late time series expansion, that describes 
the effects of inhomogeneity and anisotropy, and vanishes 
in the homogeneous isotropic FLRW limit. 

The resummed redshift and luminosity distance as a 
function of redshift are given by. 
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and, 
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where, 



f{yo,ye;k,pQ) 



ye a{yo;k,po) 
yoa{ye;k,po) 



(28) 



and where the function a{y; fc, po) may be obtained from 
the integral in ^TE\i . and we have defined, 

g{y,,Z;k,po) ^ (1 + Z) °^^°' ^' ^"^ (Z; fc, po), (29) 
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where y*{Z) is given by 
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In order to perform the resummation, we have factored 
out the quantities f{yo, ye', k, pq) and g{yo, Z; k, po). The 
advantage of using this form of the redshift and luminos- 
ity distance is that the factor involving the late time ex- 
pansion is trivial in the homogenous and isotropic limit. 
In order to compute the part which does not vanish in this 
limit we have only to determine the FLRW scale factor 
function a(yo;fc, pg) and perform the integral ([5(11). Fo'" 
practical purposes, where we wish to determine the model 
data through fitting to observations, we may simply nu- 
merically evaluate these. This should be contrasted with 
the late time expansion about flat sliced de Sitter ([24| . 
where instead one must calculate terms in the late time 
expansion even to describe the homogeneous isotropic 
limit, and truncating this expansion would lead to an 
unnecessary approximation. 

It is an interesting question whether all of the char- 
acterizing data can in principle be determined through 
observation. For example, truncating the late time ex- 
pansion to the order we have worked with above, it is 
easy to show that all components of Rij, hij and S7(iRjk) 



can be determined by extracting the quantities R, h 
and v^ffRi-j, v^v^hij and ( w™-D"w'Vi^,„„ -I- ^w*Vi^ 

from comparison with luminosity distance-redshift ob- 
servations for suitably many observation directions v^. 
However, at this order we cannot determine the non- 
symmetric parts of ViRjk, such as \/uRj-\j.. Whether 
these non-symmetric parts can be determined from 
higher orders in the expansion is a question we defer for 
later investigation. 

To conclude, by measuring the dependence of the lumi- 
nosity distance on redshift and angle on the sky (giving 
v^) an observer comoving with the matter can extract 
information about the characterizing data that lives at 
the future conformal boundary. For example, working 
to first non-trivial order they can in principle determine 
all components of the boundary Ricci tensor (and hence 
Riemann tensor) at the boundary location that they will 
asymptote to. At next order they determine all of hij at 
this location and some information about derivatives of 
the Riemann tensor there. With sufficiently good data 
going to higher orders allows an observer to extract in- 
formation about higher derivatives of gij and hij . 



V. AVERAGING SOLUTIONS 

We have used a dust fluid to model the CDM, which is 
expected to be of a particulate nature. Hence a notion of 
averaging has been introduced to arrive at a smooth fluid 
description. It is then natural to consider how the in- 
homogeneous solutions transform under averaging. Our 
geometric characterization of a solution in terms of a 3- 
geometry and its perturbation suggest that a canonical 
averaging scheme can be found by considering an aver- 
aging of 3-geometries. Ricci flow provides precisely the 
analogy of diffusion for geometry, and as such gives an 
averaging of geometry. Introducing a flow time t, we 
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have the Ricci flow equation, 



^-g^Ar) 



-2R^Mr)], 



(32) 



where we take (jij (0) to be our original metric (jij . Flow- 
ing for some time r sets the scale of smoothing. 

Now consider a metric gij and its perturbation gij + 
eSgij. Then we may Ricci flow both, and the difference 
between them after a time t will be governed by a per- 
turbation tSgijij) where, initially 5gij{Q) — Sgtj and its 
flow is given by the linearization of the Ricci flow about 

5ij(T), 



—Sgij = ALSgij - '2SIi^iUj) 



(33) 



where A^ is the Lichnerowiz operator defined by, 

Ai<55„ = V^%, + 2R/] "Sg„u^ ~ '2R^r&~9,)m , (34) 



and Ui is defined earlier in (J14p . The flows are invari- 
ant under the scaling 5gij — >• XSgij for A S M+. Thus 
we take simultaneously flowing (p2|) and ([33]) to pro- 
vide a canonical smoothing of the data characterizing 
our late time cosmology, {gij ,[5 gij]). Starting with ini- 
tial data (gjj(O), (5^^ (0)) one should perform the simul- 
taneous flows for some time T, yielding smoothed data 
{gij{T),dgij{T)). One should then perform a coordinate 
transformation to present this data in harmonic gauge 
Ui = 0, and then use this to give the solution data 
{(jij , hij ) . We emphasize that as Ricci flow is a geometric 
flow, one may start with data (^^ (0), Sgij{0)) in any co- 
ordinate system and it will always yield the same result 
for {gij, hij) up to diffcomorphisms. We note that the av- 
eraging scheme we naturally arrive at appears somewhat 
similar to the Ricci flow averaging considered by Buchert 
andCarfora [lllI3- 

Let us consider smoothing data close to a flat FLRW 
cosmology. For the flat FLRW solution gij = 6ij and 
6gij = p6ij for a constant p. This is a fixed point of our 
flow equations above. Moreover it is a stable fixed point. 
Firstly flat space is locally stable under Ricci flow, the 
flow existing for infinite time with small deformations 
'diffusing' away. Secondly consider a general perturba- 
tion Sfjij — p 6ij + Eij now for non-constant p and trace- 
less Eij, such that we have chosen harmonic gauge Ui = Q, 
so V^Eij = ^\7ip, and take p — >■ constant and Eij — ^ 
as |a;| — >■ c». 

Then about fiat space gtj = 5ij, the flow (|33p in har- 
monic gauge takes the simple form. 



d - , d - o 



(35) 



and therefore reduces to diffusion in flat space, which is 
well defined for all flow times, and renders p constant at 
late time, and Eij — )■ 0. We note that this flow preserves 
the harmonic gauge Ui = Q we started in. The physical 
constraint ft, < reduces to the condition p < 0, which 



is indeed preserved under diffusion. Thus for our pro- 
posed smoothing flows the flat FLRW cosmology model 
is locally a stable attractor, as we would physically wish. 
Since Ricci flow is a highly non-linear operation, one 
will not be able to flow arbitrary starting data for an 
infinite flow time, and singularities may develop along 
the flow. For example the closed FLRW model is not a 
fixed point of our smoothing flow since it requires gij to 
be the metric on a round sphere and this shrinks to zero 
size in a finite flow time. Note also that for our solutions 
to be physical we require the dust density to be positive, 
and in general under our above flows this may not remain 
positive for arbitrary lengths of flow. 



VI. CONCLUSION 

Suppose one wishes to consider spatial curvature in a 
cosmological model as is often done. As we emphasized 
in the introduction, at least in the context of inflation 
and presumably more generally, it is unnatural to intro- 
duce spatial curvature without considering inhomogene- 
ity and anisotropy. Once one introduces spatial curva- 
ture without assuming homogeneity and isotropy, one is 
typically in a regime beyond that described by cosmo- 
logical perturbation theory. For example, as we have 
discussed above, if the spatial curvature is introduced as 
a squashed 3-sphere, rather than a round 3-sphere, then 
for any appreciable squashing, cosmological perturbation 
theory will fail as the solution is not a small metric de- 
formation of FLRW. Thus if one wishes to observe or 
constrain the magnitude of spatial curvature in our uni- 
verse, one is naturally led to the more general problem of 
characterizing curvature, inhomogeneity and anisotropy 
non-linearly in the amplitude of the metric and matter 
field perturbations. It is this characterization, and con- 
struction of the associated cosmological solution that has 
been the purpose of this paper. 

Using Starobinsky's late time expansion (analogous to 
techniques developed for 'holographic renormalization' 
in AdS/CFT) together with a resummation of the ex- 
pansion, we are able to find the general inhomogeneous 
anisotropic solution to cosmological constant and dust 
fluid matter that is locally asymptotic to expanding de 
Sitter. The terms in the expansion describe a systematic 
deformation about a homogeneous isotropic expanding 
FLRW solution, and are arranged as a derivative expan- 
sion in the data that characterizes the solution. The 
solution is fully non-linear in the amplitude of the metric 
and matter field deformations from FLRW. The leading 
terms in the expansion describe deformations on hori- 
zon scales, and higher terms describe deformations on 
increasingly small scales. The 'final' data lives on the 
conformal boundary of the asymptotic de Sitter future. 
It is given by the boundary metric gij and the tensor hij 
that lives on that boundary geometry, and as we have dis- 
cussed, this pair can be elegantly thought of as defining 
a 3-metric and a perturbation of that 3- metric (up to an 
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inequality on the trace of the perturbation that must be 
satisfied to ensure positive dust density) . Truncating the 
expansion at some finite order then ahows a controUed 
approximation to the true inhomogeneous anisotropic so- 
lution. 

Recent cosmological observations indicate that we are 
resident in a universe whose matter content primarily 
consists of a cosmological constant and cold dark mat- 
ter. We have argued that ignoring radiation is justified 
provided we restrict attention to the recent matter era 
{z <C 10^). We expect that for weak inhomogeneity 
and isotropy, the series given by our late time expansion 
should converge for this recent matter era, and further- 
more, the convergence will be fast. Thus truncating the 
series at some low finite order will provide a good and 
controlled approximation to the true solution. Hence we 
hope that the solution is more than simply a formal con- 
struction, and may actually allow us to interpret observa- 
tions in light of our assumptions that the Einstein equa- 
tions hold and that the late time matter content is A and 
dust fluid, without assuming anything about the homo- 
geneity or isotropy of our universe. Ideally current and 
future structure observations, such as standard candle 
supernovae observations, could allow us to characterize 
systematically and precisely the large scale inhomogene- 
ity and anisotropy, measuring them or constraining their 
size. 

Our calculation of the luminosity distance-redshift re- 
lation to leading non-trivial orders in principle allows one 
to compare to large scale structure data and directly 
determine the Riemann tensor of the boundary metric 
and the value of the hij tensor at the position corre- 
sponding to the observer position. These give the first 
anisotropic corrections to the solution. It allows one to 
extract partial information on the derivative of that cur- 
vature, giving Vfi^j-fc), but not the antisymmetric com- 
ponents. This represents the leading inhomogeneity in 
the solution. Working to higher order in the expansion, 
one would extract more information about the future 
boundary data {gij , hij ) . It would be interesting to un- 
derstand whether there are parts of the data that are 
inaccessible by the luminosity distance-redshift measure- 
ment, and if so whether they might be determined by 
other types of observation, and precisely how much data 
can be extracted in principle [27|. 

The use of an effective perfect fluid description of the 
dark matter already implies that an averaging has been 
performed. We consider averaging in the context of our 
inhomogeneous and anisotropic late time solutions. We 
argue that since the flnal data can be described in terms 
of a 3-geometry and its perturbation, it is natural to con- 
sider geometric flows, such as Ricci flows, to smooth the 
data. We consider Ricci flow of the boundary metric 
gij which then induces an associated flow on perturba- 



tions of that metric, and hence on the tensor hi 



We 



confirm that using Ricci flow smoothing, a small pertur- 
bation away from flat sliced FLRW is smoothed towards 
the FLRW solution as we would wish. 
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Appendix A: Luminosity distance-redshift relation 

In order to construct the past light cone, we begin by 
considering null geodesies. We take a tangent to these 
N^ so that N^^Nfj^ = and use the y coordinate to 
parametrize the geodesic, so that N^ = 1, where the 
Greek indices fj, = {y,i}- Then defining A^* = li^/H we 
find 4-n' 

ay 



N'^d^n^ is given as. 



dy H ™" 



n n - gmnn g 



-n'n™ 1 (Al) 



which preserves the condition n^n^ gij = 1 along the null 
geodesic. The null curve is parametrized as a;*(y) so that 
n* = Hdx^{y)/dy. In the usual way, one then obtains a 
second order o.d.e. for x''(jj), which may be solved given 
some initial data constrained by n^n^gij ~ 1. Our late 
time solution ([S]) allows us to compute g^^-, gij and F- '^ 
in a ?; expansion, and hence we may construct a solution 
to the null geodesic also in such a y expansion. We now 
do this to the order at which the dust contributes in the 
y expansion. 

Consider a source, comoving with the dust, and emit- 
ting monochromatic radiation isotropically. We choose 
coordinates such that the source is located at {y,x^) = 
{ye,xl). Consider a null geodesic which passes through 
this emitter point and reaches the future conformal 
boundary at {y,x^) = (0,5;'). For calculational conve- 
nience we choose to write the spatial variation of the 
data gij{x) and hij{x) as a Taylor expansion about the 
emitter's comoving coordinate position in the following 
coordinate system. 



9i]i^) "^k 
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:R, 



kji 



(^'=-4)(^'-4) 



-1 y ■ n,-HlJ7n 

D x'-—x 



^{VuRuj^n) . Xx'^ - Xt){x' - X',){X"^ - X^) 



+ 0{X~ XeT 



h^j{x) =hi 



+ 0{x~Xe). 



This coordinate system will simplify the calculations that 
follow, although we emphasize that the results we will 
obtain are coordinate invariant. The coordinates can be 
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thought of as an extension of Riemann normal coordi- 
nates to higher order, and is discussed in an analogous 
context (null geodesies in late time expansions for pure 
cosmological constant matter) in ;2(J]. We recall that in 
3-dimensions the Riemann tensor is simply given by the 
Ricci tensor as ([TT|) . 

Later we will re-express this data in terms of the data 
at the observer's spatial position on the future conformal 
boundary. We may obtain the corresponding expansion 
for a-- to second order in x^ — x\, 






Rij 
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With data specified at x* — x\ the solution will not only 
be an expansion in y but also in y^. We will consider 
that y and j/e are of the same order in the expansion, so 
0{y) ^ 0{ye)- We employ the following ansatz. 



v' + ^-;^_{A'+yeB') + y^C^+0{y\y\ -^ 
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The data parametrizing this geodesic is w*, and is the spa- 
tial tangent to the geodesic at y = on the future confor- 
mal boundary. We make this choice of parametrization 
for ease of calculation; later we will re-express this data 
in terms of the spatial tangent at the observer position. 
Integrating this expression we obtain the comoving coor- 
dinate along the geodesic, 
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Hx^ (y) = Hxl + {y- y,y + ^L—lk (A» + y^B^) 
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where the constant of integration has been fixed by the 
condition x^{ye) = x\. With the expansions for gij,hij 
and n' the geodesic equation can be solved as an expan- 
sion in y and ye, 
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The condition gijii^n^ = 1 implies the data, v'^, is nor- 
malized as 1 = v^v^gij{y = 0) = v^v^gij{x). This yields 
the condition v^vHij = 1-1- 0{y'^) to the order required 
for our calculation. The data gij and hij is specified at 
(2/,a;^) = (0,x^). 

The redshift, Z, of a photon emitted at ye and observed 
at j/o rnay be calculated as, 
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Z = — e^ ■ 

Vo 



(A6) 



where n* solves the null geodesic equation above. By 
considering (jAll) we deduce. 
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where F"" = v'^i 

We now consider the luminosity distance. Consider 
the congruence of null geodesies emitted from the point 
y — ye, x'' = Xe, generated by taking the directions v"" — 
v^ + a5v\ -\- P5v\, where 5v\ ^ are infinitesimal vectors 
and < a,/? < 1. Ev\ 2 must be linearly independent 
of each other and v^ , and must obey Vi8v\ 2 = 0. A 
comoving dust observer sees this congruence trace out 
an infinitesimal 2-area bA, calculated as, 
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where 6x\ 2 give the variations in the position of the null 
geodesic at time y due to the direction deformations 8v\ 2, 
keeping the source position, x*, fixed. The factor of y 
preceding n* accounts for the fact that n' is normalized 
with respect to g^ rather than the spatial metric, ^. 

We now calculate the variation in a;'(y) due to varia- 
tions of the geodesic data, v^, keeping x\ fixed, 
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In order to compute (jA8[) we make use of the identity, 

e,,fe (t16v16v^ + T^Sv^SvC) = ep,nnSv^Sv^ {SftrT - Tf) , 

(AlO) 
furthermore we note that 



^5v^5vi = VpSa 
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(A5) 



for some constant Sa. To the order at which dust enters 
the expansion we find. 
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Integrating this infinitesimal area over all possible data for the null geodesic v'^ whilst keeping y fixed defines the 
area of a closed 2-surface surrounding the emitter position, A{y). Then the solid angle of the congruence in the 
rest frame of the emitter can be found from 6^ — 47rlimj,^j,^ "aTT' Then the luminosity distance, Dl, seen by 

an observer at y = y^, is defined by, D\ = (1 + Z)"^ /q° ■ In practise we compute A{y) using angular variables, 
V = cos(0)9i + sin(0) cos((/))92 + sin(6') sin((/))93 and choose 5v\ — dgv^ and (5u| = d^v^, integrating over the full range 
of angles < 6* < tt, < < 27r. We find. 
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(A13) 



or by inverting (jA7[) . we can eliminate ye and obtain luminosity distance as a function of the redshift, 

^2 72 / on , 7\ Z{l + Z) 
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The expressions for redshift (|A7[) and luminosity distance (JA14I) are characterized by data at the emitter position 
and geodesic data w*, the spatial tangent to the geodesic at the boundary, y — 0. We wish to convert these expressions 
so that all data is at the observer position. We calculate, evaluated on the geodesic at the observer position. 
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whilst data for the geodesic at (y,x') = (0,a;*) can be converted to the normalized direction vector at the observer 
position, 

v'^v' + 0{y^) (A16) 

normalized at the observer position, v^ffgij[Q) = 1. For the redshift we obtain. 
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x^=0 



where V 



mn ^ y^i-y" — ij^mn^o). For the lumiuosity distance we find. 
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